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Abstract 

In mSUGRA model we assume that gravitino, the LSP, plays the 
role of cold dark matter in the universe, while the lightest stau, the 
NLSP, catalyzes primordial BBN reconciling the discrepancy between 
theory and observations. We have taken into account all gravitino 
production mechanisms, namely decay from heavy scalar fields, de- 
cay from the NLSP, and from the thermal bath. We find that the 
dark matter constraint is incompatible with the lower bound on the 
reheating temperature. 



There is accumulated evidence both from astrophysics and cosmology 
that about 1/4 of the energy budget of the universe consists of so called 
dark matter, namely a component which is non-relativistic and does not feel 
the electromagnetic nor the strong interaction. For a review on dark matter 
see e.g. |T]. Although the list of possible dark matter candidates is long, 
it is fair to say that the most popular dark matter particle is the LSP in 
supersymmetric models with R-parity conservation [2]. The superpartners 
that have the right properties for playing the role of cold dark matter in the 
universe are the axino, the gravitino and the lightest neutralino. By far the 
most discussed case in the literature is the case of the neutralino (see the 
classical review [3]), probably because of the prospects of possible detection. 

On the other hand, primordial Big-Bang nucleosynthesis (BBN) is one 
of the cornerstones of modern cosmology. In the old days, BBN together 
with Bubble's law and CMB supported and strengthened the Hot Big-Bang 
idea. Nowadays, BBN can be used to test and constrain possible new physics 
beyond the standard model. Although the general agreement between the 
theoretical predictions and the observed light nuclei abundances is quite im- 
pressive, it is true that that there is a discrepancy between the standard 
theory and Lithium isotopes observations [1] . Recently it was proposed that 
the so called catalyzed BBN can reconcile the aforementioned discrepancy [5] . 
Consider a heavy unstable negatively charged particle X~ which can form 
Coulomb bound states together with Helium 4. Then the following reaction 



can affect the primordial light element abundances substantially. The Lithium 
6 observations can constrain the properties of X~ , and in particular its life- 
time 16]. 



In CMSSM the lightest neutralino or the lightest stau is typically the 
lightest superpartner. Therefore in scenarios in which gravitino is assumed 
to be the LSP, then the neutralino or the stau is the NSLP and therefore un- 
stable with a lifetime which is typically larger than BBN time tsBN ~ 1 sec. 
Energetic particles produced by the NLSP decay may dissociate the back- 
ground nuclei and significantly affect the primordial abundances of the light 
elements. If such processes occur with sizable rates, the predictions of the 
standard BBN scenario are altered and the success of the primordial nucle- 
osynthesis is spoiled. BBN constraints on cosmological scenarios with exotic 
long-lived particles predicted by physics beyond the Standard Model have 



{^HeX-) +D Li + X 



(1) 
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been studied [7], and the neutralino NLSP scenario is already disfavored [S]. 
However the stau NLSP is still a viable scenario, and in this case the stau 
can play the role of the heavy unstable negatively charged particle X~. 

In the present article we work in the mSUGRA model and assume that 
the gravitino is the LSP while the stau is the NLSP. Gravitino plays the 
role of cold dark matter in the universe while the stau plays the role of X~ 
to catalyze BBN. In the CMSSM, contrary to the MSSM, there is a small 
controllable number of parameters, namely four parameters and a sign. These 
are 

• Universal gaugino masses 

Ml (Mgut) = M2 (Mgut) = M3 (Mgut) = ^1/2 (2) 

• Universal scalar masses 

mj^^McuT) = mo (3) 

• Universal trilinear couplings 

A:^{Mgut) = 4(Mgc;t) = A\^{Mgut) = (4) 

• 

tanP = — (5) 

where vi , V2 are the vevs of the Higgs doublets and Mgut ~ 10^^ GeV 
is the Grand Unification scale. 

plus the sign of the n parameter from the Higgs sector. In mSUGRA there 
is the additional condition for the gravitino mass, 771.3/2 = friQ [9]. The neu- 
tralino and stau masses depend on the universal scalar and gaugino mass as 
follows P [ID] 

rriy. ~ 0.42mi/2 (6) 
m? ~ ml + 0.15ml/2 (7) 

From the above formulas we can easily see the two limits in which either 
neutralino or stau is the lightest of the usual superpartners. In the limit in 
which mo ^ ^1/2 ? the stau mass rrif ~ rriQ and in this case the neutralino 
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is lighter than stau. In the hmit in which mo <C mi/2, the stau mass ~ 
0.387mi/2 and in this case the stau is hghter than the neutrahno. The natural 
range for values of mi/2 is from 100 GeV up to a few TeV, and following |6] 
we shall consider for mi/2 the range 

100 GeV < mi/2 < 6 TeV (8) 

For the gravitino abundance we take all possible production mechanisms into 
account and impose the cold dark matter constraint |TT] 

0.075 < n^dmh^ = ^3/2h'^ < 0.126 (9) 

At this point it is convenient to define the gravitino yield, = 
where is the gravitino number density, s = heff{T)^T^ is the entropy 
density for a relativistic thermal bath, and /ig// counts the relativistic degrees 
of freedom. The gravitino abundance Q3/2 in terms of the gravitino yield is 
given by 

where we have used the values 



To = 2.73K = 2.35 x 10"^^ GeV 



(11) 

/ie//(To) = 3.91 (12) 
Pcr/h^ = 8.1 X 10"^^ GeV^ (13) 

The total gravitino yield has three contributions, namely one from the ther- 
mal bath, one from the NLSP decay, and one more from heavy moduli decay. 

y3/2 = Y,J/^ + F3%^^^ + ^372°'^"'"^ ( 14) 

The contribution from the thermal production has been computed in fl2[ \T3\ 
[13]. In [12] the gravitino production was computed in leading order in the 
gauge coupling g^, in [T3| the thermal rate was computed in leading order in 
all Standard Model gauge couplings gY,g2,g3, and in [13] new effects were 
taken into account, namely: a) gravitino production via gluon — > gluino 
+ gravitino and other decays, apart from the previously considered 2^2 
gauge scatterings, b) the effect of the top Yukawa coupling, and c) a proper 
treatment of the reheating process. Here we shall use the result of [T^ since 
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the corrections of [131 [S] do not alter our conclusions. Therefore the thermal 
gravitino production is given by 



The second contribution to the gravitino abundance comes from the decay 
of the NLSP 

' rriNLSP 

with the gravitino mass, rriNLSP the mass of the NLSP (here the stau) 
and VL^LSph'^ the NLSP abundance had it did not decay into the gravitino, 
which for the stau is estimated to be |15| 



n^h^^i^^S" (17) 
\2TeVJ ^ ' 

The decay width of stau to tau and gravitino is given by [21 [S] 



2 ^ 4 



i = r(f^r + V.,/.) = j^^ 1-^1 (18) 

where Mp = 2.4 x 10^^ GeV is the reduced Planck mass. In the limit in 
which m3/2 ^ the stau lifetime is simplified as follows 

Tf = 6.1x 10^ sec ( 4rT7 ) 19) 

\100GeVJ \ rrif ) ^ ' 

According to the results of [6], the catalyzed BBN imposes on the stau lifetime 
the upper bound 

< 5 X 10^ sec (20) 

Using that = rrio, rrif = 0.387mi/2, and the expression for the stau 

lifetime we obtain 

mi/2> 426 (^)'^' GeV (21) 

On the other hand, the catalyzed BBN only works if the lifetime of the X~ 
particle is larger than a certain value 

Tf > 10^ sec (22) 
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or 

mi/2< 588 (^)'^' GeV (23) 

For example, if m^i2 = thq = 100 GeV, then mi/2 takes values within the 
interval 

2.69 TeV < mi/2 < 3.71 TeV (24) 

The last contribution to the gravitino abundance comes from a heavy mod- 
ulus decay. Moduli are four- dimensional scalar fields with a potential, and 
they appear from higher dimensional supergravity/superstring theories upon 
compactification down to four dimensions. Moduli parameterize the shape 
and size of the manifold used for compactificationand, and typically their 
vaccum expectation value is of the order of the Planck mass Mp. Finally, 
generically and without fine-tuning their mass is close to the gravitino mass. 
For example in the KKLT model \W\ one finds that rriT — 47r^m3/2, where T 
is the total volume modulus. However, with fine tuning it is possible to have 
a heavy modulus and an arbitrarily light gravitino [17]. The gravitino yield 
from modulus decay is given by |18] 



Ymodulus _ 3/2 -L R f25) 

where Tr is the reheating temperature at the modulus decay, mx is the 
modulus mass, Vtot is the total decay rate of the modulus, and r3/2 is the 
modulus decay rate to a pair of gravitinos. Using the supergravity lagrangian 
one can compute both Vtot and r3/2- The modulus dominantly decays into 
gauge bosons and gauginos with a total decay width [T8] 



Q 3 

Ttot = r(X ^ all) ^ T{X ^ gg) + T{X ^ ^^) = — ^ (26) 

167r 

while the decay width into a pair of gravitinos is given by 

= 288^J^ ^^^^ 

in the limit in which mx ^ ^3/2- Therefore the branching ratio for the 
decay channel into a gravitino pair is given by 

Br{X ^ 7/>3/2V'3/2) = ^ = ^- 0.01 (28) 

J- tot 04 
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We assume that the modulus dominates the energy density when it decays 
and thus it reheats the universe. The reheating temperature at the modulus 
decay is determined as usual by the condition H{Tji) = Ttot, where H{T) is 
the Hubble parameter as a function of the temperature during the radiation 
era. Thus one finds 

/ on \i/4 

or 

Notice that Tr is determined entirely by the modulus mass. From BBN there 
is a lower bound on the reheating temperature, coming from the requirement 
that at BBN time all three neutrino species are thermalized |19]. This limit 
reads 

Tr>1 MeV (31) 
which in turn induces a lower bound on modulus mass 

mx > 1.5 X 10^ GeV (32) 

For given values of gravitino mass and stau mass, the total gravitino abun- 
dance is a function of the modulus mass. In figure 1 we show the grav- 
itino abundance as a function of modulus mass for mo = 100 GeV and 
mi/2 = 2.69 TeV, while in figure 2 we show the gravitino abundance versus 
the modulus mass for mo = 100 GeV and mi/2 = 3.7 TeV. In both figures the 
horizontal strip shows the allowed range for cold dark matter. The figures 
show that gravitino can explain the cold dark matter in the universe provided 
that the modulus mass is of the order of GeV, which is obviously in contradic- 
tion to the limit mx > 1.5 x 10^ GeV. Apart from that, since the modulus 
field decays into two gravitinos, kinematics requires that mx > 200 GeV. 
Therefore, for these two reasons we conclude that the scenario under inves- 
tigation must be excluded. Finally, if we consider larger values of m3/2 we 
find that the situation becomes even worse. The contribution from the NLSP 
decay increases and therefore the modulus mass must be lower than before. 

We remark in passing that if the gravitino mass is treated as one extra 
parameter and can be low enough, then it possible to satisfy both the cold 
dark matter constraint and the lower bound on the modulus mass. For 
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example, in figure 3 we show the gravitino abundance versus the modulus 
mass for a light gravitino, m^i2 = 0.1 GeV, and mi/2 = 234 GeV. We see that 
the gravitino abundance is within the allowed observational range provided 
that the modulus mass is mx ~ 10® GeV, which satisfies the limit mx > 
1.5 X 10^ GeV. This implies that if the scenario under investigation is to 
be realized in nature, then the gravitino must be light or, put it differently, 
the correct supersymmetry breaking pattern might be the gauge mediated 
one [20]. 

In summary, in the present work we have studied supersymmetric dark 
matter in mSUGRA assuming that the gravitino is the LSP and that the 
stau is the NLSP. The gravitino plays the role of cold dark matter in the uni- 
verse, while the stau catalyzes the standard BBN reconciling the discrepancy 
between theoretical predictions and observations on cosmic lithium isotopes 
observations. We have taken into account all three gravitino production 
mechanisms and have imposed the cold dark matter constraint. The pro- 
duction mechanisms for the gravitino are i) a heavy modulus decay, ii) the 
NLSP (stau) decay, and iii) scattering processes from the thermal bath. We 
have assumed that the modulus dominates the energy density when it decays 
and reheats the universe. The reheating temperature is determined entirely 
by the modulus mass. Constraints from BBN impose a lower bound on the 
reheating temperature and in turn a lower bound on the modulus mass, 
> 1.5 X 10^ GeV. On the other hand, if we compute the total gravitino 
abundance and impose the cold dark matter constraint we find that the mod- 
ulus mass should be much lower than the above lower limit. Therefore, we 
conclude that the scenario under investigation must be excluded. 
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Figure 1: Gravitino abundance versus modulus mass for mo = 100 GeV and 
mi/2 = 2.69 TeV. The strip shows the cold dark matter constraint. 
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Figure 2: Same as figure 1, but for mi/2 = 3.7 TeV. 
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Figure 3: Same as figure 1, but for 777.3/2 = 0.1 GeV and 7771/2 = 234 GeV. 
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